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ABSTRACT 


Investigations  into  methods  of  utilizing  ambient  gas 
as  propellant  at  satellite  altitudes  have  been  experimentally  and 
analytically  investigated.  As  part  of  this  investigation,  a 
class  of  space  charge  limited  plasma  accelerators  has  been  studied. 
Experimentally,  discharges  in  these  devices  are  low  current 
(.1  amp  to  2 amps)  and  high  voltage  ( 200v  to  10,000v)  . Ions 
are  produced  by  electron  bombardment  of  the  adsorbed  gases  on 
the  anode.  Treating  the  flow  of  ions  and  electrons  as  a conducting 
fluid  permits  analogies  with  fluid  flow  phenomena.  Sub-  and 
supercritical  flow  for  both  ions  and  electrons  correspond  to  sub- 
sonic and  supersonic  flow  in  a fluid.  In  one-dimensional  flow, 
the  electric  field  creates  an  "electrostatic  nozzle"  with  critical 
flow  occurring  where  the  electric  field  is  zero  ( the  "electro- 
static throat").  A rigorous  analysis  of  the  flow  gives  quali- 
tative agreement  with  measured  discharge  characteristics.  Space 
charge  buildup  limits  the  ratio  of  beam  power  to  electric  power 
input  to  under  1%  for  the  operating  mode  wherein  ions  are 
produced  from  the  adsorbed  gas  on  the  anode.  Conventional 
performance  as  a pulsed  plasma  accelerator  can  be  obtained  by 
desorbing  the  gas  from  the  anode  rapidly  so  that  breakdown  occurs. 
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1.0  INTRODUCTION 

The  thrust  to  conduct  the  maneuvers  associated  with 
satellite  reorientation,  orbit  modification  and  drag  make-up 
has  up  to  the  present  been  provided  ky  any  one  of  a number  of 
means : 

cold  gas  jets 
hydrazine  jets 
ion  engines 

arc  jets  and  plasma  thrusters  carrying 
their  own  propellant 

chemical  rocket  propulsion  systems 

colloid  thrusters. 

The  problems  and  disadvantages  of  these  systems  are  many 
and  diverse,  but  they  have  in  common  the  critical  problem  that 
the  total  impulse  available  is  proportional  to  the  mass  of 
propellant  carried  on  the  space  craft.  Mission  life-times  must 
hence  be  programmed  from  the  point  of  view  of  this  limitation, 
rather  than  from  the  mission  requirements.  Utilization  of  ambient 
atmosphere  as  propellant  eliminates  this  problem. 

Studies  have  been  conducted  to  determine  the  feasibility 
of  collecting  air  at  high  altitudes  for  use  in  chemical  or 
electrical  propulsion  systems.  The  most  extensive  investigations 
were  those  of  S.  T.  Demetriades  and  colleagues  on  the  Air-Scooping 
Orbital  Rocket  (A-SCOR)  (Ref.  1-4) . Similar  studies  were  carried 
out  elsehwere  (Ref.  5-9) . These  studies  indicated  that  the  system 
was  feasible  provided  that  a power  source  of  considerable  size 
(several  Megawatts)  was  available  to  assist  in  the  collection 
process  and  to  heat  the  air  (thermally  or  electrically)  used  as 
propellant  in  the  thruster  system.  The  lack  of  an  adequate 
power  source  appears  to  be  the  main  obstacle  to  the  imple- 
mentation of  this  system. 
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The  ideal  propulsion  system  for  satellites  would  be  one 
for  which  both  power  and  propellant  were  available  in  situ. 

Rapid  improvements  in  solar  cell  technology  and  in  the 
technology  of  axisymmetric  plasma  accelerators  over  the  past 
ten  years  make  it  desirable  at  this  time  to  investigate  the 
feasibility  of  such  a truly  infinite  total  impulse  propulsion 
system.  Obviously,  unless  extremely  large  solar  cell  arrays 
are  utilized  (on  the  order  of  square  kilometers) , the  power 
level  available  will  be  low  (under  10  kilowatts) , so  that 
the  cryogenic  collection  and  storage  of  propellant  - as  suggested 
in  the  above-mentioned  studies  - will  not  be  feasible. 

An  alternate  means  of  propellant  handling  could  result  through 
the  use  of  a space-electric  ram  jet.  The  postulated  space- 
electric  ram  jet  would  utilize  an  axisymmetric  electric 
discharge  surrounded  by  a solenoidal  magnetic  field  to 
ionize  and/or  accelerate  the  ambient  air.  When  the  power  is 
supplied  by  solar  cells  this  type  of  propulsion  system  has 
the  potential  for  extremely  large  total  impulses. 

The  space-electric  ram  jet  can  be  postulated  to 
operate  in  a number  of  modes.  The  mode  employed  can  depend 
upon  many  factors,  such  as  power  availability,  thruster 
performance  characteristics,  possible  communications 
interference,  etc.  One  possible  mode  of  operation  that  would 
appear  to  be  very  efficient  would  be  operation  of  the  thruster 
at  altitudes  where  the  ambient  ionization  level  is  near  its 
peak,  i.e.  in  the  region  between  200  and  300  km  altitude. 

By  properly  designing  the  magnetic  field  configuration,  the 
discharge  could  perhaps  encompass  a large  cross-section  in 
space  and  utilize  only  the  existing  ions  and  electrons.  In 
this  manner  the  thrust  would  be  produced  without  any 
accompanying  ionization  energy  loss.  Also,  in  this  mode,  it 
might  be  possible  to  operate  with  a high  electric  potential 
across  the  electrodes  and  hence  obtain  high  specific  impulses 
(over  3000  sec) . 
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Unfortunately,  laboratory  experiments  to  investigate 
if  such  an  operational  mode  is  feasible  are  far  beyond  the 
scope  of  the  present  investigation.  In  fact,  it  is 
likely  that  the  only  practical  test  of  this  hypothetical 
mode  of  operation  would  be  a test  in  near  space  on  an 
orbiting  laboratory,  such  as  that  envisioned  for  the  space 
shuttle . 

A steady- state  operating  mode  for  a space-electric 
ram  jet  has  been  observed  in  laboratory  experiments  at 
pressure  levels  of  over  10-3  Torr.  (Ref.  10-18)  . It  remains 

to  be  seen  if  this  mode  of  operation,  which  requires  volume 

. . . -7 

ionization,  will  perform  well  at  pressure  levels  of  10 

-3 

to  10  Torr.  Part  of  the  current  may  be  carried  in  a 
"filament"  (which  can  be  many  centimeters  in  diameter)  that 
extends  out  from  the  cathode  and  loops  back  to  the  anode. 

This  loop  encompasses  all  of  the  magnetic  flux  lines  inside 
the  anode  ring.  The  "loop"  current  rotates  rapidly 
(200-800khz)  and,  because  of  the  high  electron  concentration 
and  temperature  within  it,  acts  as  an  ionization  front, 
ionizing  the  atoms  of  the  ambient  gas  as  it  spins  through 
them.  The  ions  that  are  produced  are  then  accelerated 
(electrostatically)  toward  the  axis  and  join  with  electrons 
from  the  cathode  jet  to  form  the  exahust  beam  that  escapes 
from  the  magnetic  field  and  produces  a net  axial  thrust  ard  a 
torque  on  the  engine.  It  is  also  possible  to  consider  that 
the  thrust  is  produced  by  j x B forces,  in  which  case  the 
axial  thrust  results  from  the  interaction  of  the  azimuthal 
current  with  the  radial  component  of  the  applied  magnetic 
field. 

The  distance  that  the  current  loop  extends  out  into 
the  space  behind  the  engine  would  probably  depend  upon  the 
ambient  pressure,  the  applied  magnetic  field  strength,  and 
the  total  current.  In  some  cases  (very  low  pressures)  it 
may  be  many  meters  long  and  several  meters  in  diameter  at 
its  widest. 
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Such  an  engine  would  produce  a torque  as  well  as  a 

thrust. 

In  order  to  have  no  torque  on  the  average,  the 
magnetic  field  would  be  reversed  periodically.  Naturally, 
use  could  be  made  of  the  torque  for  spinning  or  despinning 
the  vehicle . 

It  should  be  emphasized  that  the  discharge  modes 
postulated  above  are  very  different  from  conventional  vacuum 
arcs  (Ref.  19) . In  the  latter,  the  material  which  forms 
the  plasma  and  exhaust  beam  is  almost  invariably  eroded  from 
the  cathode  or,  less  frequently,  from  the  anode,  or  insulator 
structure.  Great  care  would  have  to  be  taken  in  the  design 
of  space-electric  ram  jet  devices  to  ensure  that  this  does 
not  occur  and  that  the  major  component  of  the  discharge 
plasma  and  exhaust  beam  consists  of  air. 

A quite  different  operational  mode  can  be  postulated. 
This  would  be  a pulsed  mode  which  utilizes  the  air  adsorbed 
on  the  surfaces  of  the  engine.  A discharge  can  be  postulated 
to  ionize  and  accelerate  this  adsorbed  air  off  the  surface 
and  shape  it  into  an  exhaust  beam.  The  surface  would  then 
be  free  to  adsorb  more  air  from  the  atmosphere  and  the  cycle 
could  be  repeated. 

Of  the  three  operational  modes  postulated  above,  the 
last  is  the  one  that  lends  itself  most  readily  to  experimental 
investigation  in  the  laboratory.  Extensive  tests  have, 
accordingly,  been  conducted  on  several  engine  configurations 
that  operated  in  this  mode. 

Results  from  these  tests  have  been  reported  in 
references  ( 19  and  20) . These  results  will  be 

examined  in  detail  in  this  report  from  two  points  of  view. 
First,  the  experimental  data  will  be  compared  to  theoretical 
performance  estimates.  Second,  the  experimental  data  will 
be  used  to  generate  thruster  performance  estimates. 
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2.0 


ENGINE  CONFIGURATION  AND 
OPERATIONAL  CHARACTERISTICS 


Detailed  descriptions  of  the  measurements  made 
on  the  various  discharge  and  accelerator  configurations 
tested  can  be  found  in  references.  For  convenience  this 
information  is  summarized  here. 

2 . 1 Assembly  1. 

A special  piece  of  test  equipment  was  designed  and 
constructed  for  investigating  the  influence  upon  the  discharge 
characteristics  of  desorption  effects  on  the  electrode 
surfaces  of  an  axisymmetric  plasma  accelerator.  A picture 
of  this  device  is  shown  in  Figure  1.  It  consists  of  2 NRC 
ionization  gauges  (type  807)  facing  each  other  and  connected 
together  by  a glass  cylinder  two  inches  in  diameter  and  seven 
inches  long.  A glass  cylinder  three  inches  in  diameter 
and  twelve  inches  long  is  connected  to  this  assembly  to 
increase  the  volume.  A Veeco  ionization  gauge  is  connected 
to  the  system  for  measuring  pressure.  This  whole  assembly 
was  mounted  in  a furnace  and  connected  to  a vacuum  pumping 
system.  The  system  was  baked  at  a temperature  of  about 

350°C  for  several  days  while  it  was  evacuated  to  a pressure 
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of  under  10  Torr  at  which  time  it  was  sealed  off.  The 
volume  of  assembly  1,  is  2.56  litres  (156  cubic  inches). 

One  of  the  NRC  tubes  was  used  an  as  axisymmetric 
plasma  accelerator.  This  was  accomplished  by  using  the 
cylindrical  collector  plate  as  an  anode  and  the  heated 
filament  as  the  cathode.  The  grid  was  allowed  to  float 
electrically.  A coil  of  300  turns  of  number  18  wire  was  wound 
around  the  outside  of  the  gauge  to  produce  a solenoidal  magnetic 
field.  This  coil  was  two  inches  long,  had  an  I.  D.  of  two 
inches  and  an  0.  D.  of  four  inches.  The  coil  was  calibrated 


and  found  to  produce  32.5  gauss  per  ampere  of  current  on  the 
center  line. 
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The  filament  was  heated  by  discharging  a 60,000  y.  f. 
condenser  bank  through  it.  The  bank  could  be  charged  to 
several  hundred  volts,  but  it  was  found  that  charging  it  to 
voltages  of  25  to  31  volts  provided  adquate  heating.  The 
scope  traces  shown  in  Figure  2 show  the  current  through  the 
filament  and  the  thermal  radiation  from  the  filament  as 
measured  by  a photo  cell,  both  as  a function  of  time.  It 
is  seen  that  the  peak  filament  temperature  is  attained  about 
60  milli-seconds  after  the  current  through  the  filament 
is  initiated.  One  side  of  the  filament  was  grounded. 

The  magnet  coil  was  energized  by  another  60,000  y.  f. 
condenser  bank.  The  scope  trace  in  Figure  3 shows  both  the 
current  through  the  magnet  coil  and  the  current  through  the 
filament  as  a function  of  time.  The  magnetic  field  strength 
has  a time  constant  for  decay  of  about  70  milli-seconds. 

Energy  for  the  discharge  between  the  ring  anode  and 
the  heated  filament  was  supplied  from  a condenser  charged  to 
a high  voltage . Tests  have  been  conducted  with  condenser 
banks  of  2,  4 and  32  y.  f.  charged  to  voltages  of  from  1,000 
volts  to  8,000  volts. 

The  other  NRC  ionization  gauge  was  used  as  an  ion 
detector.  The  circuit  shown  in  Figure  4 was  used.  The 
differential  amplifier  was  used  so  that  both  electrodes  of 
the  probe  could  be  left  floating  to  ensure  that  no  current 
flowed  between  the  discharge  tube  and  the  ion  detector. 

2.1.1.  Summary  of  Experimental  Results  Using 

Assembly  1. 

Over  500  independent  tests  were  conducted  with  this 
piece  of  equipment.  A unique  type  of  slow  discharge  has  been 
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observed  in  this  device.  Typical  current-voltage  traces 

from  an  oscilloscope  are  shown  in  Figure  5.  At  high  voltages 

— 4 

(over  8,000  volts)  and  high  ambient  pressure  (over  10  Torr) 
breakdown  occurred  and  a more  conventional  arc  discharge 
resulted.  The  peak  current  in  the  "slow"discharge  was  usually 
a fraction  of  an  ampere  and  the  discharge  lasted  tens  of 
milli-seconds . The  current  transistioned  discontinuously 
(in  less  than  2 y.  sec.)  to  a much  lower  value  as  the  potential 
drop  across  the  electrodes  reached  some  low,  critical  value 
( * 1,200  volts) . 


2 . 1 . 1 . 1 Influence  of  Ambient  Pressure 

It  was  possible  to  reduce  the  pressure  in  the  assembly 
— 8 

to  1 x 10  Torr  by  using  the  Veeco  ionization  gauge  as  an 
ion  pump  for  a period  of  several  days.  Also,  by  heating  up 

the  grid  of  the  Veeco  gauge  as  well  as  the  tube  it  was  possible 

_ 3 

to  raise  the  pressure  in  the  assembly  to  10  Torr.  Tests 

— 3 • 8 

were  hence  conducted  over  a pressure  range  of  10  to  10 
Torr. 

No  magnetic  field  was  applied.  The  filament  condenser 

bank  was  charged  to  30  volts  and  a 4 y.  f.  condenser  was 

charged  to  4,000  V for  the  discharge.  Current  and  voltage  traces 

were  recorded  as  a function  of  time  using  an  oscilloscope.  (Fig.  5) 

Negligible  changes  in  the  discharge  characteristics  occurred 

— 8 —5 

as  the  pressure  was  raised  from  3 x 10  Torr  to  5 x 10  Torr. 
From  this  pressure  on,  however,  as  the  pressure  was  increased 
the  following  trends  were  observed: 

- the  peak  current  increased  with  the  pressure. 

the  magnitude  of  the  discontinuous  current 
change  increased. 

the  voltage  drop  at  which  the  discontinuous 
current  change  occurred  decreased  as  the 
pressure  was  increased. 
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2. 1.1. 2 Influence  of  Charging  Voltage 

on  4 u.  f.  Discharge  Capacitor 

A series  of  tests  were  conducted  in  which  the  pressure 
was  held  constant  at  3.5  x 10-6  Torr . No  "discharge" 
occurred  until  the  charging  potential  was  over  1200  volts. 
Increasing  the  voltage  on  the  capacitor  beyond  this  point 
increased  both  the  peak  current  and  the  duration  of  the 
current.  Characteristic  current-voltage  relationships  are 
plotted  in  Figure 


2.1.1. 3 Influence  of  Size  of  Discharge  Capacitor 


Many  tests  were  conducted  with  a 16  u.  f.  capacitor 
charged  to  4000  volts.  In  general,  the  peak  current  was 
only  slightly  higher  than  when  a 4 y.f.  was  used.  The  time 
duration  of  the  discharge  was,  however,  5-6  times  as  long. 

In  some  cases  the  current  discontinuity  was  observed,  and  in 
others  not.  The  discharge  lasted  long  enough  so  that  the 
heating  current  through  the  filament  had  decayed  to  a small 
fraction  of  its  peak  value.  Since  at  times  the  voltage  on 
the  discharge  capacitor  did  not  fall  below  3000  volts,  it  can 
be  concluded  that  there  is  not  sufficient  ion  bombardment 
of  the  cathode  to  sustain  the  discharge. 


2. 1.1. 4 Influence  of  an  Applied  Solenoidal 

Magnetic  Field 

When  a solenoidal  magnetic  field  was  applied  to  the 
discharge  a number  of  effects  were  observed.  The  discharge 
current  decreased  as  the  magnetic  field  was  increased  and 
was  eventually  suppressed  completely  at  field  strengths  of 
several  thousand  gauss.  Oscillations  in  the  discharge  current 
occurred  at  moderate  magnetic  field  strengths  (several 
hundred  gauss)  at  a frequency  of  about  1,000  hz. 
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A very  interesting  interaction  between  the  cathode 
* emission  and  the  magnetic  field  was  observed.  The  magnet 

j coil  was  excited  simultaneously  with  the  filament  current 

and  a peak  magnetic  field  of  about  1,350  gauss  occurred 
after  a few  milli-seconds . The  energy  to  he;at  the  filament 
was  then  decreased  successively  for  each  test  point.  The 
discharge  current  rose  as  this  happened  and  the  oscillations 
died  out.  At  a quite  low  energizing  energy  for  the  filament 
the  discharge  current  was  eventually  suppressed. 

2 . 2 Assembly  2 

Because  of  the  small  volume  of  Assembly  1 a strong 
possibility  existed  that  some  interaction  could  occur  between 
the  discharge  and  the  walls  of  the  vessel.  Also,  the 
■ chamber  was  so  small  that  no  accurate  measurement  of  a time 

of  flight  could  be  made  to  obtain  some  estimate  of  the 
particle  velocity  in  the  exhaust  beam.  Accordingly,  a series 
of  tests  were  planned  and  conducted  in  a large  volume  facility, 
using  the  same  engine  configuration  as  in  the  earlier  tests 
and  covering  the  same  range  of  operation  parameters. 

] 

A vacuum  chamber  located  at  the  Convair  Aerospace 
Division  of  General  Dynamics  in  San  Diego  was  used  for  the 
tests.  The  specifications  for  the  chamber  "A"  are  given  in 
Table  I . Because  of  the  time  and  cost  involved  in  pumping 
I down  and  repressurizing  the  chamber,  it  was  decided  to 

accomplish  as  many  tests  as  possible  with  only  one  pumpdown. 
This  required  that  special  precautions  be  taken  to  ensure 
that  no  malfunction  of  any  part  of  the  electrical  or 
mechanical  equipment  in  the  tank  could  occur.  Duplicate 
engines  were  mounted  in  the  tank  and  four  ion  detectors 
were  positioned  downstream  of  the  source.  The  placement 
i of  these  components  are  shown  in  Figure  6.  A photograph 
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of  the  duplicate  source  (engine)  is  shown  in  Figure  7.  The 
ion  detectors  were  all  identical  and  a photograph  of  one  is 
shown  in  Figure  8. 
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A great  deal  of  care  was  taken  to  maintain  the  noise 
level  in  the  detector  circuits  at  a very  low  level.  All 
electrical  connections  between  the  tank  and  the  instrumentation 
were  made  with  #8  shielded  amphenol  cable.  A special  feed- 
through plate  was  designed  and  constructed  which  ensured 
that  shielding  continued  through  the  tank  walls  and  also 
that  no  ground  loops  could  occur  in  the  tank.  A photograph 
of  the  feed-thru  plate  is  shown  in  Figure  9. 

All  of  the  electrical  and  diagnostic  equipment  needed 
for  the  test  was  assembled  in  a mobile  laboratory  from  which 
all  tests  were  conducted  on  the  site.  A schematic  of  the 
electric  circuit  for  the  ion  source  (engine)  is  shown  in 
Figure  10.  All  the  shields  for  the  cables,  the  chassis  and 
the  instrumentation  were  grounded  on  a common  bus  inside  the 
mobile  laboratory  to  minimize  ground  loop  effects. 

The  circuits  for  the  ion  detectors  were  designed  to 
measure  accurately  picoamperes  of  beam  current  impinging 
on  the  detector  faces.  The  circuits  used  are  shown 
schematically  in  Figure  11.  The  detector  consisted  of  a 
thin  wall  outer  cylindrical  nickel  electrode  3 cm  in  diameter 
by  2 cm  long  and  a central  solid  copper  electrode  1 cm  in 
diameter  and  2 cm  long.  These  electrodes  were  rigidly  mounted 
in  a plexiglass  plate  and  the  electrical  connections  to  the 
leads  were  buried  in  the  plexiglass.  The  polarity  was  such 
that  the  outer  electrode  would  collect  ions  and  the  inner 
one  of  the  electrons. 

Note:  Throughout  the  testing  period  the  test  environment 
was  kept  at  liquid  nitrogen  temperature. 
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2.2.1 


Summary  of  Experimental  Results  Using  Assembly  2 


During  a period  of  3 days  over  300  data  points  were 
taken.  A discharge,  similar  to  that  observed  in  Assembly  1, 
was  observed.  Breakdown  occurred  much  more  readily, however , 
hence  most  of  the  tests  were  conducted  with  the  discharge 
capacitor  charged  to  only  2000  volts. 

2.2.1. 1 Influence  of  Ambient  Pressure 

Data  was  taken  on  the  discharge  characteristics  and 
probe  signals  as  the  tank  pressure  was  bled  up.  During  these 
tests  a standard  "accelerator"  operating  condition  was 
used,  where  a peak  solenoidal  magnetic  field  of  several 
hundred  gauss  occurred  in  the  region  of  the  electrodes. 

No  significant  change  in  the  discharge  characteristics  was 
observed  over  the  pressure  range  of  1 x 10  ® Torr  up  to 
3 x 10~5  Torr.  At  higher  pressures  high  frequency  oscillations 
were  observed  and  the  peak  current  decreased. 

These  oscillations  were  measured  on  the  detectors 
as  well  as  across  the  discharge.  The  magnitude  of  the 
discontinuous  change  in  the  current  remained  constant,  but 
the  potential  drop  at  which  the  discontinuity  occurred 
decreased  as  the  pressure  was  increased  above  3 x 10  ^ Torr. 

2. 2. 1.2  Influence  of  Changing  Voltage  on 

4 p.f.  Discharge  Capacitor 

A series  of  tests  were  conducted  at  the  minimum 
pressure  obtainable  near  the  beginning  of  the  test  series. 

The  voltage  on  the  discharge  capacitor  was  increased  in 
increments  of  500  volts,  starting  at  1000  volts.  The  peak 
current  increased  rapidly  with  the  initial  voltage  across 
the  discharge;  however,  the  time  duration  of  the  discharge 
increased  only  marginally.  Breakdown  occurred  intermittently 
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when  the  charging  potential  was  over  2500  volts.  Discharge 
characteristics  are  displayed  in  Figure  12. 

2. 2. 1.3  Influence  of  an  Applied  Solenoidal 

Magnetic  Field 

As  the  solenoidal  magnetic  field  strenth  was  increased, 
the  peak  current  of  the  discharge  decreased.  The  time 
duration  of  the  discharge  increased  with  the  magnetic  field. 

At  modestly  high  magnetic  field  strengths,  oscillations 
occurred  in  the  current  and  voltage  of  the  discharge. 

Increasing  the  magnetic  field  strength  eventually  led  to  a 
breakdown. 

2 . 3 Discussion  of  Experimental  Results 

The  insensitivity  of  the  "slow"  discharge  charact- 
eristics to  the  ambient  pressure  leads  to  the  conclusion 
that  the  discharge  is  associated  with  ions  produced  on  the 
anode  surface  due  to  electron  bombardment.  Examination 
of  the  discharge  curves  also  leads  to  the  conclusion  that 
the  processes  are  dominated  by  space  charge  effects. 

The  measured  signals  on  the  detectors  in  the  experiments 
on  Assembly  2 indicate  that  plasma  has  been  accelerated 
downstream  of  the  discharge.  Two  approaches  towards  estim- 
ating the  performance  potential  of  the  device  as  a space 
propulsion  system  become  available: 

i)  attempt  to  measure  the  thrust,  the  mass  flow  rate 
and  the  electrical  power,  or  alternatively  the 
impulse, the  mass  utilized  per  pulse  and  the 
electrical  energy.  Use  these  measured  quantities 
to  calculate  the  specific  impulse  and  the  efficiency 
of  the  device. 
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ii)  identify  the  mechanisms  and  processes  from 
the  experimental  results  and  construct  a 
theoretical  model  to  predict  performance 
capabilities . 

Past  history  of  plasma  propulsion  studies  indicates 
that  approach  (i)  can  often  lead  to  over-optimistic 
performance  estimates  due  to  the  difficulty  of  the  measure- 
ments and  to  interactions  occurring  between  the  engine  and 
the  test  environment.  Where  feasible,  then,  approach  (ii) 
has  much  to  recommend  it.  The  flow  of  non-colliding 
charged  particles  can  be  accurately  modelled  and  a well 
defined  boundary  value  problem  can  be  formulated  and 
solved.  This  approach  was  hence  taken  and  the  analysis  of 
the  charged  particle  flow  (treated  as  a continuum)  is 
presented  in  the  next  section. 
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3.0 


ANALYSIS 


The  flow  of  a completely  ionized  plasma  with  or 
without  currents  needs  to  be  investigated  in  order  to 
determine  the  magnitude  of  the  ion  and  electron  flux  rates 
in  the  projected  engine  configurations.  It  is  appropriate 
to  examine  the  flow  in  the  limit  when  no  volume  collisions 
between  the  ions  and  electrons  occur.  This  is  a problem 
that  has  been  studied  by  numerous  investigators,  starting 
with  Langmuir.  The  innovation  attempted  here  is  to  utilize 
the  fluid  flow  equations. 

In  order  to  gain  some  understanding  of  the  procedure 
and  of  the  significance  of  the  results,  the  problem  is 
investigated  in  three  phases: 

Phase  I one-dimensional  ion  or  electron 
flow  - the  classical  Langmuir- 
Childs  problem. 

Phase  II  one-dimensional  ion  and  electron 

flow.  Each  species  flow  constitutes 
a current. 

quasi-one-dimensional  flow  of  ions 
and  electrons. 


Phase  III 
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PHASE  I 


ONE  DIMENSIONAL  EQUATIONS  FOR  ELECTRON  OR  ION  FLOW 
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Mass  conservation: 

| e | n w ■ - j = constant. 


Momentum 
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conservati on : 
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Energy  conservation: 
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Poisson's  equation: 

P-  = M" 

**  eo 


Equation  of  State: 
p * n k T 
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3.1.2  \ 
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THERMODYNAMIC  AND  GASDYNAMI C RELATIONS 
FROM  MOMENTUM  AND  ENERGY  CONSERVATION 
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becomes 


3.1.6 


From  equation  of  state 


3.1.7 


From  continuity 
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WORKING  EQUATIONS 


Introduce  (8)  into  (3) 
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Introduce  (1 ) into  (4) 
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3.1.10 


Multiply  (19)  by  E and 
Introduce  (9)  into  (10) 
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Equations  (9)  and  (11),  together  with  boundary 
conditions  permit  a complete  solution  to  the 
probl em. 


0 ( Langmui r-Chi 1 ds  Solution) 


SOLUTION  FOR  hc  * 


Boundary  Conditions: 


uc  * o 
Ec  « 0 

vc  ■ 0 


Integration  of  eq.  (9)  gives: 
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Integration  of  eq.  (11)  gives 
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Combining  eq.  (12)  and  (13)  and  integrating  again 
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Rewriting  in  conventional  form 
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GENERAL  SOLUTION 

1)  Eq.  (9)  indicates  that  for  u to  increase 

monatoni cal ly  with  x the  electric  field  must 
be  positive  for  u < u and  negative  for  u > u 
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4)  Rewrite  equations  (9)  and  (11) 
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3.1.11a 


5)  Integrate  (11a)  from  the  critical  point  to  an 
arbitrary  point. 
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6) 


Introduce  E into  (9a). 
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3.1.19 


using  eq.  17  and  introducing  a non-dimensional  space 
variable  E,  an  equation  that  can  be  evaluated  in  terms 
of  the  conditions  at  the  cathode  is  obtained. 
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Once  a solution  is  obtained  by  numerical  integration,  it  is 
possible  to  compare  this  more  exact  solution  with  the 
Langmuir-Childs  solution.  This  is  done  by  deriving  an 
expression  that  is  unity  in  the  limit  of 


“Tef  <<:  ^ " Vc 


This  expression  is: 


■,<**  - Ve>  V^A-V 


9/4  £a2  Mc1/3 


MA2-Mc2  +3{^ 


L - C.  C. 


An  expression  for  the  voltage  can  be  derived  in  terms  of  the  cathode 
temperature  and  the  "Mach  Numbers"  at  the  cathode  and  anode. 

YkT  ? M y-1  Y-l  2 2 

va  ■ VC  = y^t  {( } '1}  + Mc  {ma  ‘ Mc 


in  evaluating  this  expression  a cathode  temperature  of  3000  K 
was  used. 

T * 3000° K 
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A reference  distance  can  be  computed  by  using  the  maximum 
emission  capability  of  tungsten  at  3000°K  as  the  "Mach  1" 
condition  to  evaluate  nc> 
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COMPARISON  WITH  LANGMUIR  CALCULATIONS 
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Phase  II 

No  Collisions  between  Ions  and  Electrons 
One-dimensional  flow. 
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| 3.2  Phase  II 

, No  Collisions  between  Ions  and  Electrons 

One-dimensional  flow. 
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Similarly  for  the  ions 
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The  two  momentum  equations , 3.2.2  and  3.2.5  can  be  added 
together  and  integrated  to  give: 
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Using  equations  3.2.11  and  3.2.13 
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or 
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Equations  3.2.14b  and  3.2.15a  permit  us  to  evaluate  the 

"terminal  velocities  U.**  and  U * in  terms  of  the  ratios  of 

i e 

the  electron  and  ion  pressures  and  temperatures  at  the  two 
"critical"  positions. 
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Phase  III 

Plasma  and  Current  Flow 
No  Collision  Interaction 
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Derive  critical  condition  for  ions 
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from  continuity 
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Momentum  equation  becomes 
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To  eliminate  the  ion  temperature,  use  entropy  condition  that 
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Next,  eliminate  E by  multiplying  3.3.11  and  A/A*  and 
differentiating.  First,  define 


Collecting  terms  and  using  equation  3.3.9,  we  obtain  a second 
order  differential  equation  for  the  ion  velocity 
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Equation  3.3.12  couples  the  ion  and  electron  velocities. 
Therefore,  a differential  equation  for  the  electron  velocity 
must  next  be  derived.  Proceeding  as  above  for  the  ions,  the 
electron  momentum  equation  becomes 
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The  electron  temperature  can  now  be  eliminated  by 
introducing  the  entropy  constraint  for  the  electrons,  namely 
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3.3.13 


The  electric  field  can  be  eliminated  from  this  equation  in  the 
same  manner  as  it  was  in  the  equation  for  the  ion  velocity. 

The  result  is  a second  order  differential  equation  for  the 
electron  velocity.  For  convenience,  define 
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The  differential  equation  for  the  electron  velocity  becomes 


In  some  cases,  it  is  more  convenient  to  use  a relation  between 
the  ion  and  electron  velocities 


'i 


The  a lr>ve  equations,  together  with  the  appropriate  boundary 
conditions  permit  solutions  to  a wide  variety  of  problems. 


Boundary  Conditions 


Some  boundary  conditions  are  obtained  at  the  position  where  the 
ion  velocity  is  critical.  These  are  listed  below 
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A similar  set  of  boundary  conditions  are  obtained  where  the 
electron  velocity  becomes  critical,  these  are  listed  below 
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3.4 


Procedure  to  Evaluate  Solution  for 
Counter-flowing  Ions  and  Electrons. 


Choose  values  for  the  following  quantities 


i)  The  non-dimensional  potential  drop 
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ii)  The  charge  on  the  ion,  Z. 

iii)  The  ratio  of  the  critical  values  of  the 
electron  to  ion  temperature, 
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Using  the  electron  energy  equation,  by  iteration, 
evaluate  the  non-dimensional  electron  velocity  at 
the  position  where  the  ion  velocity  is  critical. 
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3)  Using  the  ion  energy  equation,  by  iteration,  evaluate 

the  non-dimensional  ion  velocity  at  the  position 
where  the  electron  velocity  is  critical. 
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3.4.4 


4)  Make  a first  choice  of  the  value  of  vhe 

"Childs-Langmuir" , constant  L.C. 
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5)  Evaluate  the  ratio  of  the  values  of  the  ion  pressure 

and  electron  pressure  at  their  respective  critical 
positions 
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6)  Evaluate  the  charge  density  ratio 
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7)  Numerically  integrate  Poisson's  Equation, 
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Using  the  energy  rela  ions  as  follows: 

a)  when  0 > n < 1 ; U +U.  >1. 

e l 

iterate  for  Ug  and  using  equations. 
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b)  when  n<0  U.>  1;  U < 1. 
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iterate  for  IK  using  Equation  3.4.10. 


iterate  for  Ug  using 
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c)  when  n > 1/  IK<  1;  Ug  >1. 


iterate  for  Ug  using  Equation  3.4.9 


iterate  for  IK  using 
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Compare  electric  field  from  integration  with  value 
calculated  from  the  momentum  equation. 
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3.5  Relations  Between  Actual  and 
Non-Dimensional  Variables 
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Some  typical  solutions  of  the  equations  are  presented  in 
Figures  13  through  21.  A graph  of  the  "Childs-Langmuir" 
coefficient  as  a function  of  the  potential  difference 
between  the  two  critical  positions  is  presented  in  Figure  22. 
There  are  a number  of  interesting  aspects  to  the  solutions: 


1)  Attempts  to  obtain  non-symmetrical  solutions  have 

been  unsuccessful.  If,  indeed,  only  symmetrical 
solutions  exist,  then  the  critical  points  will  be 
located  where 

Z kTg**  = kTj,*  and  the  following  relations  hold 
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2)  In  the  "sub-critical"  region,  the  pressure  and 

temperature  ratios  between  the  "critical"  position 
and  the  source,  usually  a hot  surface,  can  be  very 
large . 


This  can  be  seen  in  Figures  15  and  18,  where  the 
logarithm  of  the  temperature  and  pressure  ratios 
are  plotted  as  a function  of  the  non-dimensional 
distance . 


As  expected,  the  region  of  influence  of  the  ambient 
or  surface  temperature  upon  the  flow  concentrates 
into  smaller  regions  as  the  ratio  of  the  potential 
drops  to  the  temperature  increases. 


Interelectrode  potential  drops  can  be  much  lower  than 
the  potential  drop  between  the  two  critical  positions. 

The  charge  neutralization  that  occurs  due  to  the  ion 
current  flow  permits  the  electron  current  to  be  up 
to  an  order  of  magnitude  higher  than  occurs  with  space 
charge  limited  electron  flow  only. 
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4. 


COMPARISON  OF  EXPERIMENTAL  AND  ANALYTIC  RESULTS 


The  analytic  results  qualitatively  agree  with  the  ex- 
perimental results  as  summarized  in  Fig. 22  . The  results 
cannot  be  expected  to  agree  quantitively  since  the  experiments 
were  conducted  with  equipment  that  had  cylindrical,  rather  than 
planar  flow.  Although  the  analytic  procedures  for  calculating 
cylindrical  counter-flow  of  ions  and  electrons  are  presented 
in  this  report,  no  solutions  have  been  yet  calculated.  Another, 
more  interesting  phenomenon  that  has  not  yet  been  analytically 
treated  is  the  development  of  criteria  for  transition  from  the 
ion-electron  flow  mode  to  the  electron  flow  mode  only.  According 
to  the  experiments,  this  transistion  occurs  abruptly.  It  can  be 
speculated  that  this  occurs  when  the  ion  density  at  the  anode 
drops  to  a value  where  p*  <.  p**  . The  critical  flow  for  the  ions 
cannot  then  be  established  and  the  electric  field  at  the  anode 
rises  abruptly. 

The  following  description  of  the  mechanisms  in  the  "slow" 
discharge  follows  from  the  experimental  and  analytic  studies 
described  earlier: 

1.  A "discharge"  that  is  dominated  by  space  charge 
effects  has  been  obtained.  The  characteristis  of  the 
discharge  can  be  analytically  investigated  by  study 
of  the  behaviour  of  counter-flowing,  non-colliding 
electron  and  ion  beams. 

2 . The  discharge  is  initiated  by  application  of 
several  kilovolts  between  concentric  electrodes,  and 
heating  the  thermionic  cathode,  which  is  centrally 
located . 

3 . The  space  charge  limited  electron  flow  produces 
ions  on  the  anode  surface  by  inelastic  collisions  with 
the  adsorbed  surface  layers  of  gas. 

4 . A space  charge  limited  flow  of  ions  towards  the 
cathode  occurs.  This  flow  partially  neutralizes  the 
cathode  space  charge  and  permits  a much  higher  electron 
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current  density  than  occurs  with  electron  flow  only. 

5 . Both  the  electron  and  ion  current  density  are 
controlled  by  a "critical  flow"  criterion  at  an 
electrostatic  throat,  rather  than  by  emission  or  ion 
production  limitations. 

6 If  the  power  density  at  the  electrodes  becomes 
too  high,  gas  is  desorbed  at  a high  rate.  The  gas 
density  between  the  electrodes  can  then  become  high 
enough  to  permit  breakdown. 

7 . When  the  ion  density  near  the  anode  surface  drops 
to  some  critical  value,  the  ion  flow  through  the 
"electrostatic"  throat  cannot  be  established  and 
space  charge  limited  electron  flow  only,  results. 

8 . The  ratio  of  ion  current  to  electron  current  is 
equal  to  the  square  root  of  the  ratio  of  the  electron 
to  ion  masses. 


5.0  IMPLICATIONS  OF  THE  ANALYTIC  AND  EXPERIMENTAL  RESULTS 

FOR  THRUSTER  PERFORMANCE 

Considerable  insight  into  the  performance  characteristics 
of  electromagnetic  accelerators  utilizing  ambient  gas  as  the  pro- 
pellant has  been  obtained  from  the  experimental  and  analytic 
results  obtained  in  this  program.  This  information  is  summarized 
in  the  following. 

5.1  MODES  OF  OPERATION 

5.1.1  Mode  I Acceleration  of  Ambient  Ionospheric  Ions 

This  study  could  not  address  this  mode  of  operation 
experimentally  due  to  facility  and  cost  limitations.  Analytically, 
a correct  solution  to  a properly  formulated  boundary  value  problem 
for  this  type  of  system  is  still  several  steps  away  from  an  orderly 
investigation  of  space  charge  limited  flows.  Further  development 
of  the  approach  initiated  in  Section  3.3  of  this  report  should  per- 
mit accurate  prediction  of  performance  potential.  If  this  looks 
promising,  then  experiment  on  the  space  shuttle  would  probably  be 
the  most  cost  effective  method  of  evaluating  actual  performance. 

5.1.2  Mode  II  Acceleration  of  Ions  Produced  from  Adsorbed 
Gas  Layers  on  or  Near  the  Anode  Surface . 

The  experimental  and  analytic  investigations  of  this  program 
have  been  primarily  associated  with  this  mode  of  operation,  as 
exempified  by  the  "slow  discharge"  described  earlier.  In  terms  of 
conventional  performance  parameters,  this  mode  of  operation 
results  in  very  inefficient  use  of  the  electrical  power.  The  high 
electron  current  relative  to  the  ion  current  results  in  an 
"efficiency"  of  under  1%,  where  it  is  defined  as  beam  power  or 
energy  divided  by  electrical  input  power  or  energy.  The  specific 
impulse  of  the  exhaust  beam  in  this  mode  of  operation  can  be 
arbitrarily  high.  The  ions  are  accelerated  through  most  of  the 
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potential  drop  between  the  electrodes,  which  has  been  experiment- 
ally varied  between  ten  thousand  and  several  hundred  volts. 

5.1.3  Mode  III  Acceleration  of  Ions  Produced  by  Volume 
Collisions  of  Gas  Desorbed  from  Engine  Structure. 

This  mode  of  operation  is  equivalent  to  that  of  conventional 
pulsed  plasma  accelerators,  and  has  been  extensively  studied  by 
many  experimental  and  analytic  groups  (Refs.  25-30  ) using 

injected  or  ablated  propellant.  Performance  will  be  comparable 
to  that  observed  in  such  studies,  and  hence  no  detailed  investiga- 
tions of  this  mode  were  made. 

5.1.4  Mode  IV  Acceleration  of  Ions  Produced  by  Volume 
Collisional  Interaction  of  an  Extended  Discharge  Into 
the  Ambient  Gas. 

This  mode  of  operation  had  earlier  been  observed  (Ref.  31  ) 

in  many  electromagnetic  accelerators  with  ambient  pressures  above 

10  ^ Torr.  It  was  hoped  that  operation  in  this  mode  could  be 

obtained  at  much  lower  pressures  by  having  the  discharge  encompass 

very  large  volumes  of  gas . The  tests  conducted  in  the  large  tank 

at  General  Dynamics/San  Diego  indicated  that  volume  ionization 

effects  start  to  influence  the  discharge  only  a pressures  of 
-3  -4 

at  least  10  Torr  to  10  Torr,  approximately.  (Initially,  as 
pressure  is  increased  above  these  levels,  phenomena  occur  which  are 
probably  associated  with  a moving  or  oscillating  ionization  front.) 

As  pressures  over  10~^  Torr  occur  at  altitudes  below  60  miles  and 
hence  are  of  minimal  interest  for  satellite  propulsion, 

-3 

studies  were  not  conducted  at  pressures  above  10  Torr,  and  no 
further  information  or  data  generated  in  this  study  on  this  operating 
mode . 
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Photo-cell  .1  volts/cm 

Filament  current  200  amps /cm 

Sweep  speed  20  m.  sec. /cm 

V ■ 30  volts 


Fig.  2 Radiation  from  Filament  & Current  Through 
the  Filament. 
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Figure  4 Ion  Detector  Circuit  Using  NRC507  Ionization  Gauge. 
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Figure  9 Photogrnoh  of  Feed-Through  Plate  into  Vacuum  Chambe 
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Date  1-29-75 


Experiment 


Picture  # J 


Capacitor  Charging  Voltages 
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Lower  Detector  V 
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Upper  Trace 

Signal  Discharge  Voltage 
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Figure  12a.  Characteristic  Current  Voltage  Traces  — 
Applied  Magnetic  Field. 
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Picture  # 64  A 
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Picture  # 64  b 
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Figure  I2g  (Continued) 
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FIGURE  13 


SCHEMATIC  OF  GEOMETRIC  AND  ELECTRIC  CONFIGURATION 
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ELECTRIC  FIELD  DISTRIBUTIONS 


FIGURE  17  NON-DIMENSIONAL  ELECTRON  AND 
ION  VELOCITY  DISTRIBUTIONS 


FIGURE  18  NON-DIMENSIONAL  TEMPERATURE 
AND  PRESSURE  DISTRIBUTIONS 


FIGURE  19  NON-DIMENSIONAL  POTENTIAL  AND 
ELECTRIC  FIELD  DISTRIBUTIONS 


FIGURE  20  NON-DIMENSIONAL  ELECTRON  AND 
ION  VELOCITY  DISTRIBUTIONS 


log(p  /p**)  log(p./p*) 


FIGURE  22  NON-DIMENSIONAL  POTENTIAL  AND 
ELECTRIC  FIELD  DISTRIBUTIONS 


FIGURE  23  NON-DIMENSIONAL  ELECTRON  AND 
ION  VELOCITY  DISTRIBUTIONS 


ELECTRON  CURRENT  DENS I' 


DISCHARGE  CURRENT  (AMPERES  ) 

FIGURE26  DISCHARGE  CHARACTERISTICS,  NO  B-FIELD 


GENERAL  DYNAMICS 

Convair  Aerospace  Division 


TABLE  I 

Chamber  B was  used 


DESCRIPTION 

SPACE  SIMULATOR,  SPACECRAFT  TEST  LABORATORY 

General  Dynamics  Convair  Aerospace  Division 
San  Diego,  California 


Internal  Dimensions 


Chamber  "A"  12  ft.  diam,  19  ft.  long  horizontal  cylinder. 
Chamber  "B"  11  ft.  diam,  22  ft.  long  horizontal  cylinder. 


Pumps 

Mechanical 

Diffusion 


Traps 

Net  Speed  (from  chamber) 


3 to  6,  280  CFM  each 


Chamber  "A” 
Chamber  "B" 


Chamber  "A" 
Chamber  "B” 
Chamber  "A" 
Chamber  "B” 


3-15,000  1/s,  mercury  diffusion 
2-95,000  1/s,  oil  diffusion 
1-9500  CFM,  oil  diffusion  (booster) 
-40°F  freon  plus  LN2  Chevron 
chilled  H2O  manifold  plus  LN  Chevron 
6000  1/s  2 

80,000  to  95,000  1/s 


Chamber  Air  Conditioning 
Vacuum 


Data 


Location 


500  CFM,  -10°F  dewpoint,  60-70°F,  absolute  filters. 

Chamber  "A"  1 x 10“^  Torr  ultimate  pressure 
Chamber  "B"  <1  x 10“^  Torr  ultimate  pressure 
•with  cold  wall 

A 200  guarded-channel  data  system  is  capable  of  recording 
on  either  or  both  punched  tape  and  printed  tape  at  5 points 
per  second.  Computers  and  plotters  are  available  for  data 
reduction.  Guarding  and  isolation  of  instrumentation  pro- 
vides high  accuracy  accruisitlon  of  millivolt  signals  from 
thermocouples,  etc.  100  guarded  copper  constantan  and  an 
additional  100  copper  wire  passthroughs  are  available. 

Thirty  land  lines  to  the  Analog  and  Digital  Computer  Labors* 
tories  provide  system  control  or  additional  data  handling. 

General  Dynamics  Convair  Aerospace,  Spacecraft  Test 
Laboratory,  Bldg.  28,  San  Diego,  California,  Telephone 
277-8900.  Ext.  2221. 


